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Elliptic flow of the hot, dense system which has been cre- 
ated in nucleus-nucleus collisions develops as a response to 
the initial azimuthal asymmetry of the reaction region. Here 
it is suggested that the magnitude of this response shows a 
"kinky" dependence on the centrality of collisions for which 
the system passes through a first-order or rapid transition 
between quark-gluon plasma and hadronic matter. We have 
studied the system Pb(158AGeV) on Pb employing a recent 
version of the transport theoretical approach RQMD and find 
the conjecture confirmed. The novel phase transition signa- 
ture may be observable in present and forthcoming experi- 
ments at CERN-SPS and at RHIC, the BNL collider. 



One of the most important goals of the heavy-ion pro- 
grams in the ultrarelativistic energy domain is the search 
for the phase transition between hadronic matter and the 
quark-gluon plasma (QGP). After a decade-long effort 
based on numerous experiments at fixed-target machines 
(CERN-SPS, BNL-AGS) heavy-ion physics can be con- 
sidered a mature field today. Thus it may seem surpris- 
ing that there still is a shortage of reliable signatures 
for the elusive state QGP and the transition itself |Q. 
It may be easier to comprehend the difficulties to iden- 
tify the QGP if one takes into account that properties 
of the QGP and the transition have to be reconstructed 
from the final state which obviously is of hadronic na- 
ture. On a more fundamental level, it has become clear 
only recently that the properties of strongly interacting 
matter even far above the critical temperature T c are 
essentially non-perturbative. This makes many of the 
"first-generation" QGP signals which are based on per- 
turbation theory unreliable at best. 

On the other side, information about the QGP and 
the phase transition region has become available with 
the advent of more powerful lattice gauge simulations of 
quantum chromodynamics (QCD) ^||. Most notably, it 
has been shown that chiral symmetry is restored at rather 
low temperatures (in the range 140 to 170 MeV). Further- 
more, the Equation of State (EOS) varies rather rapidly 
in the transition region. It is not clear yet whether the 
transition is of weak first-order or just a rapid cross-over 
between the two phases. The EOS extracted from the lat- 
tice clearly displays the transition from hadron to quark- 
gluon degrees of freedoms. Pressure and energy density 
approach the ideal Stefan-Boltzmann values at tempera- 
tures > 3T C . A generic feature of the EOS in the transi- 



tion region is the presence of the so-called "softest point 
of the EOS" j|,H related to the effect that the energy 
density may jump with increasing temperature but the 
pressure does not. 

The collective transverse flow which develops in the 
heavy-ion collisions reflects on the properties of the EOS. 
Usually, one distinguishes various types of transverse 
flow, the radial (isotropic component), directed (sideward 
kick in the reaction plane) and the elliptic flow, the lat- 
ter being a preferential emission either along the impact 
parameter axis or out of the reaction plane (squeeze-out) 
H . The general idea why a phase transition may show up 
in flow observables is rather straightforward. At densi- 
ties around the softest point the force driving the matter 
expansion gets weakened. A long time ago, van Hove 
has suggested that the multiplicity dependence of aver- 
age transverse momenta may display a plateau and a 
second rise [Q. So far, it has not been possible to deduce 
the presence of a phase transition from the transverse 
momentum spectra. Some time ago we have suggested 
that the elliptic flow may be a better-suited observable 
to identify a first-order type phase transition [|). Here 
we make good on this promise and present a novel signa- 
ture of the QCD phase transition. We predict a rather 
characteristic centrality dependence of the elliptic flow if 
the created system passes through the softest region of 
the EOS in the heavy-ion reactions. 

Elliptic flow in the central region of ultrarelativistic 
collisions is driven by the almond-shape of the participant 
matter in the transverse plane ||. It was argued in || 
that elliptic flow may be more sensitive to the early pres- 
sure than the isotropic radial flow. "Early" and "late" is 
defined by the time scale set by the initial transverse size 
r t = \J (x 2 + y 2 ) of the reaction region. Thus early flow 
appears at times « r t /c while we would refer to flow gen- 
erated at times > 2r t /c as late. One reason for the larger 
sensitivity of the elliptic flow to early pressure is that the 
generated flow asymmetry works against its cause and 
diminishes the spatial asymmetry on a time scale propor- 
tional to \/ (y 2 ) — v (x 2 ). Furthermore, the elliptic asym- 
metry is proportional to the difference between the flow 
strengths in x (parallel to impact parameter) and y direc- 
tion. Thus it is more fragile than radial flow. Viscosity 
related non-ideal effects tend to wash out the pressure- 
driven asymmetries. Obviously, these effects will be more 
pronounced in the later dilute stages of the reaction. Un- 
fortunately, this could not be demonstrated in the earlier 
work. The transport model RQMD (version 2.3) |[o| em- 
ployed for the calculations lacked any sizable transverse 
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pressure in the early stages - a combination of softness 
from pre-equilibrium motion and absence of a QGP phase 
which would generate more pressure than the resonance 
matter simulated in the model. As a result the final 
hadron momentum spectra showed azimuthal asymme- 
tries much smaller than hydrodynamical results which in- 
clude a phase transition into the QGP. In the mean time, 
NA49 has analysed data for semi-central Pb(158AGeV) 
on Pb collisions jll| . The measured azimuthal asymme- 
tries are roughly equally distant from the closest results 
based on hydrodynamics and from the RQMD calcula- 
tions [|iJ,[k|. Both of these calculations show a factor 
of two disagreement, however, in different directions. In 
this Letter we are going to present results from calcula- 
tions with a new version of the transport model RQMD 
(version 3.0) which incorporates an EOS with 1st order 
phase transition. Comparing to the results obtained in 
the model without QGP phase we may assess the impor- 
tance of the phase transition. 



15 



10 



CL 

\ 
CD 




nance gas 



quark— gluon plasma 



0.5 1 1.5 2 2.5 3 3.5 4 
e (GeV/fm 3 ) 

FIG. 1. Equation of states implemented into RQMD: ra- 
tio of energy density e divided by pressure p. The dashed 
line represents the resonance gas EOS, the solid line the EOS 
including a first order phase transition with T c = 160 MeV. 

Let us first shortly describe how the phase transition 
is implemented into the model. A detailed description 
of the algorithm will be presented elsewhere. In RQMD 
nucleus-nucleus collisions are calculated in a Monte Carlo 
type fashion. While the nucleons from each colliding 
nuclei pass through each other, they are decomposed 
into constituent quarks. Strings may be excited, and 
overlapping strings fuse into ropes (with larger chromo- 
electric field strength). After their decay and fragmen- 
tation secondaries emerge and may interact with each 
other. Formed resonances are treated as unstable quasi- 
particles. This leads to a resonance gas EOS if there are 
no corrections from other interactions. The QCD dynam- 
ics in the phase transition region is not well understood. 
Even if there is a quasi-particle description it is not ob- 
vious which one of the many possible choices (strings, 
constituent quarks, partons, either massless or with dy- 



namical masses) is to be prefered. In this situation we 
have decided to stick to the implemented degrees of free- 
dom and modify the collision term instead. Since we 
expect that hydrodynamics is a reasonable approach for 
the transverse dynamics in the ultradense stage, the EOS 
should be the most relevant ingredient for the expansion 
dynamics anyway. It is well-known that different treat- 
ment of interactions between quasi-particles may modify 
the EOS. In general, if particles are free between interac- 
tion points the virial theorem specifies that the pressure 
of the system in equilibrium is given by p4[ 



P = Pi, 
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The first term arises from free streaming. The second 
term represents the non-ideal contribution AP due to 
momentum changes Sp a at discrete collision points r a . 
d=3 is the number of space dimensions, V the volume 
of the system, AT a sufficently large time interval, and 
the summation goes over all collisions, a and b speci- 
fies which quasi-particles collide. The standard collision 
term in RQMD is manufactured such that it does not 
contribute to the pressure. Now, we depart from this 
"ideal" collision term and let each quasi-particle inter- 
act elastically with a neighbor after any of the standard 
collisions. The momentum change is constrained by 

(Sp a r a ) + (Sp b n) = d ■ — ■ (AC + At s b c ) . (0.2) 
p 



At^ c refers to the time which has elapsed since the last of 
the EOS modulating collisions, p is the density of quasi- 
particles. Introducing collisions according to eq. (0.2) 
changes the pressure of the system to Pid + AP. Eq. (0.2) 
provides a numerically rather efficient method to modify 
the ideal EOS. The physics content of eq. d0.2j ) is that 
the momentum transfer may be chosen to be either suit- 
ably repulsive (QGP at high temperature) or attractive 
(mixed phase). Fig. 1 displays the ideal EOS based on 
counting the propagating quasi-particles in RQMD. In 
addition, an EOS is shown which may be produced by 
introducing energy den sity dependent additional interac- 
tions according to eq. (0.2). This EOS is the one which 
will be used for the calculations presented in this Letter. 
It is calculated from a quasi-particle model of quarks and 
gluons with dynamical thermal masses |l^,[l{|. We have 
chosen this EOS, because it provides a good fit to lat- 
tice data. The EOS contains a 1st order transition at 
T c =160MeV with a latent heat of 467 MeV/fm 3 . For 
the RQMD calculations of nucleus-nucleus collisions the 
novel interaction term is introduced in a lo cal density ap- 
proximation, i.e. all variables in eq. ( p.2[ ) are evaluated 
in the local rest system of the energy current. Neither 
is the modulation of the local pressure tensor restricted 
to regions of local equilibrium nor is - the other extreme 
- any local equilibration enforced, e.g. by randomizing 
directions of local momenta. 
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FIG. 2. Time evolution of transverse momentum 

anisotropy parameter a calculated with RQMD for the sys- 
tem Pb(158AGeV) on Pb at two impact parameters. Only 
quasi-particles within a central rapidity cut (±0.7) are in- 
cluded. 

Let us now turn to results of RQMD calculations which 
contain a phase transition. We have chosen to do the 
calculation for the system Pb(158AGeV) on Pb, i.e. the 
heavy-ion reaction at highest beam energy which is cur- 
rently available. This may be a good place to look for 
the phase transition. The time evolution of the azimuthal 
asymmetry parameter a (momentum flow asymmetry) 

«=(<p2>-<pS>)/(<i£> + <pS» (o-3) 

for quasi-particles around midrapidity is displayed in 
fig. 2. It shows a very different behaviour than corre- 
sponding calculations based on RQMD without QGP- 
type EOS B|. We see from the figure that in case of 
QGP formation essentially all of the finally observable 
asymmetry is created at times smaller than 4 fm/c. The 
mixed phase leads to a marked dip of the asymmetry for 
more central collisions. Since the pressure is comparably 
low, free motion between interactions is able to destroy 
some of the earlier created flow asymmetry. On the other 
side, the calculations for semi-peripheral collisions (e.g. 
b=7.6 fm) show that the softening in the mixed phase 
cannot stall the expansion of the system. Needless to 
say that this is a function of the latent heat which is 
very moderate for the chosen EOS. The overall effect of 
mixed phase and purely hadronic stage is rather small in 
a broad impact parameter range. Under the condition of 
an already well-developed flow asymmetry, diffusive pro- 
cesses and thermal pressure driven work seem to neutral- 
ize each other at the later stages. In the QGP-scenario 
the azimuthal asymmetry is indeed mostly a signature of 
the early pressure. It is amuzing that non-ideal effects 
from viscosity in the low-density stage may be helpful to 
infer information about the pressure in the high-density 
region. 

In the following we will present the main result of 
the Letter, the measurable azimuthal asymmetry of fi- 
nal hadrons which the experimentalists usually take to 
be the number flow asymmetry v 2 

v 2 = (cos(2$) (0.4) 



as a function of centrality. Tight impact parameter cuts 
can be obtained using a forward-energy trigger like NA49 
has. Of course, the spatial asymmetry of the reaction 
zone which is correlated with the asymmetry of the par- 
ticipant nucleons in the ingoing nuclei 

a x ={(y 2 )-(x 2 ))/((x 2 ) + (y 2 )) . (0.5) 

is itself a function of the impact parameter. Trivially, %)% 
goes to zero for very small and very large impact parame- 
ters. The value for v 2 at any given centrality reflects both 
the strength of the spatial asymmetry and the response of 
the created system due to the generated pressure. How- 
ever, we may disentangle the effects from geometry and 
dynamics. In general, the final flow asymmetry v 2 can be 
viewed as a function of many variables, a x , the average 
initial energy density eo, the transverse size r t , to name 
just a few: 

V2 = f(a x ,e ,r t , ...) » A 2 (a^) ■ a x + 0((a x -a^) 2 ) , 

(0.6) 

where we have obtained the second equation from a Tay- 
lor expansion around some intermediate value 7577, and 
taking into account that v 2 vanishes for a x — > 0. In 
Pb on Pb collisions, a x varies between and 0.50 for 
impact parameters less than 12 fm. Picking an interme- 
diate value of 7x7; means that the neglected higher order 
terms in (a x — a^) are expected to be rather small in 
practice, on the order of 10 percent. Defining the scaled 
flow asymmetry as 

A 2 = v 2 /a x (0.7) 

will therefore allow to assess the dynamical response of 
the created system to the initial spatial asymmetry. 

We display the scaled flow asymmetry A 2 versus im- 
pact parameter b in fig. 3. Of course, the asymmetry fac- 
tor A 2 will tend to vanish in the most peripheral collisions 
(6 w 2Rpb). Initial energy densities are too small, and 
the system size does not sustain extended reaction times. 
Both for pions and for protons, A 2 shows a pronounced 
variation for smaller b values. This is a result of the 
EOS softness at intermediate energy densities. However, 
non-equilibrium effects, in particular partial thcrmaliza- 
tion initially and system-size dependent freeze-out, also 
play a major role. Extracted A 2 values from hydrody- 
namic calculations || show essentially no centrality de- 
pendence, except for the grazing collisions. This feature 
is in marked contrast to the transport calculation which 
includes the non-equilibrium aspects of the dynamics. 

Without phase transition, the asymmetry factor A 2 
calculated from RQMD would simply increase monoto- 
neously with centrality - approximately linearly with the 
initial system size in the reaction plane (~ 2Rpt, — b). In- 
deed, the hard QGP stage of the reaction leads to a rapid 
increase of the asymmetry A 2 in collisions with b >10 fm 
as is visible from fig. 3. In this range of centralities the 
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initial source size \J (x 2 ) along the impact parameter axis 
is small enough that the associated characteristic time for 
the development of flow falls within the life time of the 
QGP phase. In somewhat more central collisions further 
increase of the asymmetry is cut-off after the system en- 
ters into the stage of soft and lateron viscous expansion. 
Initial energy densities change less with increasing cen- 
trality than the system size. Therefore at the character- 
istic time for flow development typical energy densities 
are in the region of the softest point. In these reactions 
increasing reaction time which is helpful to develop the 
asymmetry is counteracted by the softness of the mat- 
ter. In any case, the centrality dependence of the flow 
asymmetry follows a different slope than for the class of 
more peripheral collisions. For collisions with b < 5 fm 
kinetic equilibration which takes place on a scale of 3-4 
fm/c may be realized already in the QGP phase. This 
gives rise to yet another centrality dependence of the flow 
asymmetry A 2 (a second rise). 



the influence of a 1st order transition on elliptic flow is 
also being discussed. 
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FIG. 3. Scaled azimuthal asymmetry parameter 

A2=V2/a x of protons and pions as a function of impact pa- 
rameter for the same system and acceptance window as in 
fig. 2. 

van Hove's original idea to look for a plateau and 
a second rise in momentum spectra as a signal of the 
QCD phase transition may turn out to be true, after all. 
Present experience tells that it will probably not been 
found in the multiplicity dependence of average trans- 
verse momenta. He did not take into account that the 
dynamics of the hadronic stages may add a late radial 
flow component Jl7j which spoils the original idea. How- 
ever, if the presented calculations contain some truth it is 
much better justified to neglect the late hadronic stages 
for the azimuthal asymmetries of particle spectra. The 
presented calculation contains some uncertainties. The 
equation of state is not well-determined in the transition 
region. The admixture of baryons in the central region 
and the strong pre-equilibrium deformation of the local 
stress tensor add to the uncertainties. Nevertheless, the 
potential rewards in terms of insight into the phase tran- 
sition dynamics should justify a careful search for struc- 
ture in the centrality dependence of elliptic flow at SPS 
and future RHIC energy. 
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